Search for: All records

Determining the fate of subducted oceanic crust is critical for understanding material cycling through Earth’s deep interior and sources of mantle heterogeneity. A key control on the distribution of subducted slabs over long timescales is the bridgmanite to post-perovskite phase transition in the lowermost mantle, thought to cause rheological weakening. Using high-resolution computational models, we show that the ubiquitous presence of weak post-perovskite at the core-mantle boundary can facilitate or prevent the accumulation of basaltic oceanic crust, depending on the amount of weakening and the crustal thickness. Moderately weak post-perovskite ( ~ 10–100× weaker) facilitates segregation of crust from subducted slabs, increasing basalt accumulation in dense piles. Conversely, very weak post-perovskite (more than 100× weaker) promotes vigorous plumes that entrain more crustal material, decreasing basalt accumulation. Our results reconcile the contradicting conclusions of previous studies and provide insights into the accumulation of subducted crust in the lowermost mantle throughout Earth’s history.

 

The amount of vapor in the impact-generated protolunar disk carries implications for the dynamics, devolatilization, and moderately volatile element isotope fractionation during lunar formation. The equation of state (EoS) used in simulations of the giant impact is required to calculate the vapor mass fraction (VMF) of the modeled protolunar disk. Recently, a new version of M-ANEOS (Stewart M-ANEOS) was released with an improved treatment of heat capacity and expanded experimental Hugoniot. Here, we compare this new M-ANEOS version with a previous version (N-SPH M-ANEOS) and assess the resulting differences in smoothed particle hydrodynamics (SPH) simulations. We find that Stewart M-ANEOS results in cooler disks with smaller values of VMF and in differences in disk mass that are dependent on the initial impact angle. We also assess the implications of the minimum “cutoff” density (ρ_c), similar to a maximum smoothing length, that is set as a fast-computing alternative to an iteratively calculated smoothing length. We find that the low particle resolution of the disk typically results in >40% of disk particles falling toρ_c, influencing the dynamical evolution and VMF of the disk. Our results show that the choice of EoS,ρ_c, and particle resolution can cause the VMF and disk mass to vary by tens of percent. Moreover, small values ofρ_cproduce disks that are prone to numerical instability and artificial shocks. We recommend that future giant impact SPH studies review smoothing methods and ensure the thermodynamic stability of the disk over simulated time.

 

Abstract Large igneous provinces (LIPs) have been linked to both surface and deep mantle processes. During the formation, tenure, and breakup of the supercontinent Pangea, there is an increase in emplacement events for both continental and oceanic LIPs. There is currently no clear consensus on the origin of LIPs, but a hypothesis relates their formation to crustal emplacement of hot plume material originating in the deep mantle. The interaction of subducted slabs with the lowermost mantle thermal boundary and subsequent return-flow is a key control on such plume generation. This mechanism has been explored for LIPs below the interior of a supercontinent (i.e., continental LIPs). However, a number of LIPs formed exterior to Pangea (e.g., Ontong Java Plateau), with no consensus on their formation mechanism. Here, we consider the dynamics of supercontinent processes as predicted by numerical models of mantle convection, and analyse whether circum-supercontinent subduction could generate both interior (continental) and exterior (oceanic) deep-mantle plumes. Our numerical models show that subduction related to the supercontinent cycle can reproduce the location and timing of the Ontong Java Plateau, Caribbean LIP, and potentially the Shatsky Rise, by linking the origin of these LIPs to the return-flow that generated deep mantle exterior plumes. 

Abstract. Due to the increasing availability of high-performance computing over the past few decades, numerical models have become an important tool for research in geodynamics.Several generations of mantle convection software have been developed, but due to their differing methods and increasing complexity it is important to evaluate the accuracy of each new model generation to ensure published geodynamic research is reliable and reproducible.Here we explore the accuracy of the open-source, finite-element codes ASPECT and CitcomS as a function of mesh spacing using low to moderate-Rayleigh-number models in steady-state thermal convection.ASPECT (Advanced Solver for Problems in Earth's ConvecTion) is a new-generation mantle convection code that enables modeling global mantle convection with realistic parameters and complicated physical processes using adaptive mesh refinement (Kronbichler et al., 2012; Heister et al., 2017).We compare the ASPECT results with calculations from the finite-element code CitcomS (Zhong et al., 2000; Tan et al., 2006; Zhong et al., 2008), which has a long history of use in the geodynamics community.We find that the globally averaged quantities, i.e., root-mean-square (rms) velocity, mean temperature, and Nusselt number at the top and bottom of the shell, agree to within 1 % (and often much better) for calculations with sufficient mesh resolution.We also show that there is excellent agreement of the time evolution of both the rms velocity and the Nusselt numbers between the two codes for otherwise identical parameters.Based on our results, we are optimistic that similar agreement would be achieved for calculations performed at the convective vigor expected for Earth, Venus, and Mars. 

Phase transitions play an important role for the style of mantle convection. While observations and theory agree that a substantial fraction of subducted slabs and rising plumes can move through the whole mantle at present day conditions, this behaviour may have been different throughout Earth’s history. Higher temperatures, such as in the early Earth, cause different phase transitions to be dominant, and also reduce mantle viscosity, favouring a more layered style of convection induced by phase transitions. A period of layered mantle convection in Earth’s past would have significant implications for the secular evolution of the mantle temperature and the mixing of mantle heterogeneities. The transition from layered to whole mantle convection could lead to a period of mantle avalanches associated with a dramatic increase in magmatic activity. Consequently, it is important to accurately model the influence of phase transitions on mantle convection. However, existing numerical methods generally preclude modelling phase transitions that are only present in a particular range of pressures, temperatures or compositions, and they impose an artificial lower limit on the thickness of phase transitions. To overcome these limitations, we have developed a new numerical method that solves the energy equation for entropy instead of temperature. This technique allows for robust coupling between thermodynamic and geodynamic models and makes it possible to model realistically sharp phase transitions with a wide range of properties and dynamic effects on mantle processes. We demonstrate the utility of our method by applying it in regional and global convection models, investigating the effect of individual phase transitions in the Earth’s mantle with regard to their potential for layering flow. We find that the thickness of the phase transition has a bigger influence on the style of convection than previously thought: with all other parameters being the same, a thin phase transition can induce fully layered convection where a broad phase transition would lead to whole-mantle convection. Our application of the method to convection in the early Earth illustrates that endothermic phase transitions may have induced layering for higher mantle temperatures in the Earth’s past.

 

Mantle convection models based on geophysical constraints have provided us with a basic understanding of the forces driving and resisting plate motions on Earth. However, existing studies computing the balance of underlying forces are contradicting, and the impact of plate boundary geometry on surface deformation remains unknown. We address these issues by developing global instantaneous 3‐D mantle convection models with a heterogeneous density and viscosity distribution and weak plate boundaries prescribed using different geometries. We find that the plate boundary geometry of the Global Earthquake Model (GEM, Pagani et al., 2018,https://doi.org/10.1177/8755293020931866), featuring open plate boundaries with discrete lithospheric‐depth weak zones in the oceans and distributed crustal faults within continents, achieves the best fit to the observed GPS data with a directional correlation of 95.1% and a global point‐wise velocity residual of 1.87 cm/year. A good fit also requires plate boundaries being 3 to 4 orders of magnitude weaker than the surrounding lithosphere and low asthenospheric viscosities between 5 × 10¹⁷and 5 × 10¹⁸ Pa s. Models without asthenospheric and lower mantle heterogeneities retain on average 30% and 70% of the plate speeds, respectively. Our results show that Earth's plate boundaries are not uniform and better described by more discrete plate boundaries within the oceans and distributed faults within continents. Furthermore, they emphasize the impact of plate boundary geometry on the direction and speed of plate motions and reaffirm the importance of slab pull in the uppermost mantle as a major plate driving force.

 

Abstract This paper provides an overview of the new features of the finite element library deal.II, version 9.4. 

Abstract This paper provides an overview of the new features of the finite element library deal.II, version 9.3. 

SUMMARY Mantle convection and long-term lithosphere dynamics in the Earth and other planets can be treated as the slow deformation of a highly viscous fluid, and as such can be described using the compressible Navier–Stokes equations. Since on Earth-sized planets the influence of compressibility is not a dominant effect, density deviations from a reference profile are at most on the order of a few percent and using the full governing equations poses numerical challenges, most modelling studies have simplified the governing equations. Common approximations assume a temporally constant, but depth-dependent reference profile for the density (the anelastic liquid approximation), or drop compressibility altogether and use a constant reference density (the Boussinesq approximation). In most previous studies of mantle convection and crustal dynamics, one can assume that the error introduced by these approximations was small compared to the errors that resulted from poorly constrained material behaviour and limited numerical accuracy. However, as model parametrizations have become more realistic, and model resolution has improved, this may no longer be the case and the error due to using simplified conservation equations might no longer be negligible: while such approximations may be reasonable for models of mantle plumes or slabs traversing the whole mantle, they may be unsatisfactory for layered materials experiencing phase transitions or materials undergoing significant heating or cooling. For example, at boundary layers or close to dynamically changing density gradients, the error arising from the use of the aforementioned compressibility approximations can be the dominant error source, and common approximations may fail to capture the physical behaviour of interest. In this paper, we discuss new formulations of the continuity equation that include dynamic density variations due to temperature, pressure and composition without using a reference profile for the density. We quantify the improvement in accuracy relative to existing formulations in a number of benchmark models and evaluate for which practical applications these effects are important. Finally, we consider numerical aspects of the new formulations. We implement and test these formulations in the freely available community software aspect, and use this code for our numerical experiments.